Vacuum energy induced by an impenetrable flux tube 

of finite radius 



o 



X 



Volodymyr M. Gorkavenko * 

Department of Physics, Taras Shevchenko National University of Kyiv, 
6 Academician Glushkov ave., Kyiv 03680, Ukraine 



S ; Yurii A. Sitenkoj Olexander B. Stepanov* 



Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine, 



O . 14-b Metrologichna str., Kyiv 03680, Ukraine 



Q_i' Abstract 

d : 

. We consider the effect of the magnetic field background in the form of a tube of 

the finite transverse size on the vacuum of the quantized charged massive scalar field 
which is subject to the Dirichlet boundary condition at the edge of the tube. The 
vacuum energy is induced, being periodic in the value of the magnetic flux enclosed 
in the tube. Our previous study in J. Phys. A: 43, 175401 (2010) is extended to 
the case of smaller radius of the tube and larger distances from it. The dependence 
of the vacuum energy density on the distance from the tube and on the coupling to 
the space-time curvature scalar is comprehensively analyzed. 



1 Introduction 

The energy which is induced in the vacuum of quantized matter fields that are subject 
to boundary conditions has been studied intensively over more than six decades since 
Casimir [1] predicted a force between grounded metal plates, see reviews in Refs. [2] and 
[3J. The induced vacuum energy in bounded spaces gives rise to a macroscopic force 
between bounding surfaces. The Casimir force between grounded metal plates has now 
been measured quite accurately and agrees with his predictions, see, e.g. Refs. [I] and [5], 
as well as other publications cited in Refs. [2] and [3j. 

In the present paper we study the vacuum energy which is induced by boundary 
conditions in space that is not bounded but, instead, is not simply connected, being an 
exterior to a straight infinitely long tube. This setup is inspired by the famous Aharonov- 
Bohm effect [6], and we are interested in polarization of the vacuum which is due to 
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imposing a boundary condition at the edge of the tube carrying magnetic flux lines inside 
itself. 

Throughout the present paper, we restrict ourselves to the case of quantized scalar 
matter. A peculiarity of this case is that the energy-momentum tensor depends on the 
coupling (£) of the scalar field to the scalar curvature of space-time even then when space- 
time is flat. If scalar field is massless, then conformal invariance of the theory is achieved 
at f = £ c , where [7] 



and d is the spatial dimension; note that £ c varies from to 1/4 when d varies from 1 to 
oo. Up to now the study was restricted to the case of a singular magnetic vortex only [10] 
[13], i.e. when the transverse size of the flux-carrying tube is neglected. Therefore, the 
aim of the present paper is to take account of nonzero transverse size of the flux-carrying 
tube (for a preliminary study, see Ref. [H]). 



2 Vacuum energy density 

The operator of the quantized charged scalar field is represented in the form 




A 



where a\ and a\ {b\ and 6a) are the scalar particle (antiparticle) creation and destruction 
operators satisfying commutation relations; wave functions ^a( x ) form a complete set of 
solutions to the stationary Klein-Gordon equation 

(_V 2 + m 2 )^A(x) = ^(x), (3) 

V is the covariant derivative in an external (background) field and m is the mass of 
the scalar particle; A is the set of parameters (quantum numbers) specifying the state; 
E\ = E-\ > is the energy of the state; symbol J) denotes summation over discrete and 

A 

integration (with a certain measure) over continuous values of A. 

We are considering the static background in the form of the cylindrically symmetric 
magnetic vortex of finite thickness, hence the covariant derivative is V = d — ieV with 
the vector potential possessing only one nonvanishing component given by 

V v = $/2tt (4) 

outside the vortex; here $ is the vortex flux and ip is the angle in the polar (r, ip) coordi- 
nates on a plane which is transverse to the vortex. The Dirichlet boundary condition on 
the edge (r = vq) of the vortex is imposed on the scalar field: 

^U o =0, (5) 

i.e. quantum matter is assumed to be perfectly reflected from the thence impenetrable 
vortex. Provided the orthonormalization condition is satisfied, 

J d 3 x^x' = (A|A'), (6) 
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the solution to and (JSJ) in the case of the impenetrable magnetic vortex of thickness 
2r takes form 

i> knkz (x) = (2ny 1 e ik ^e in ^ n (kr )x 

X [Y\n-e>i>/2n\{kro)J\ n - e $/2n\{kr) - J| n -e$/27r| (^o)^jn-e0/27r| (kr)} , (7) 

where z is the coordinate along the vortex, 

f3 n (kr Q ) = [Y^ /M {kr Q ) + J^ /27rl (kr )]- 1/2 , (8) 

and < A; < oo, — oo < fc 2 < oo, n 6 Z (Z is the set of integer numbers); J^{u) and Y^{u) 
are the Bessel functions of order /i of the first and second kinds. 

In general, the vacuum energy density is determined as the vacuum expectation value 
of the time-time component of the energy-momentum tensor, that is given formally by 
expression 

£ = ( vac | [d ^d ^ + d Vd & - (£ - 1/4) V 2 (^ + tftft)] |vac) = 

= ^a^(x) ^a(x) - (£ - 1/4) V 2 VI (x) ^a(x). (9) 

A A 

In the following we shall restrict our consideration to the plane z = which is orthogonal 
to the vortex. 

Thus, the renormalized vacuum energy density in the case of the finite-thickness vortex 
takes form 

{oo 
J dkk (k 2 + m 2 ) 1/2 [S(kr, kr ) - S(kr, fcr )|$=o] - 
o 
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-(£- 1/4) A / dkk(k 2 + m 2 ) 1/2 [S(kr,kr )- S(kr,kr )\$ =0 }\ , (10) 

where, in view of ([7|), 
S(kr,kr ) = J2P 2 n(kr )x 



\Y\n-e§/2n\{kra)J\ n - e <i,/2-K\{kr) - J^-e^^ikr^Y^^^ikr)] 2 , (11) 



and A = d 2 + r~ 1 d r is the transverse radial part of the laplacian. 

Owing to the infinite range of summation, the last expression is periodic in flux $ 
with a period equal to 2ne~ 1 , i.e. it depends on quantity 



2tt 



2^ 



(12) 



where [[«]] is the integer part of quantity u (i.e. the integer which is less than or equal to 
«)■ 
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Let us rewrite (fTTl) in the form 

S(kr, kr ) = S (kr) + Si(kr, Ax ), 



(13) 



where So(kr) corresponds to the appropriate series in the case of the vacuum polarization 
by a singular magnetic vortex [11] - [T3] : 



kr 



S (kr) = £ [Jn +F (kr) + J 2 +1 - F (kr)] = dr [J F (r)J_ 1+F (r) + J-f(t)Ji-f(t)] , (14) 

n=0 { 

and Si (Ax, Axo) is a correction term due to the finite thickness of a vortex: 

oo 

Si (At, kr ) =2J2 



n=0 



t n v v v J„+F(A;ro)y;+F(A;r) - r n+F (Ax )J„+F(Ax) 
J n+ F(/cr )r n+ F(A;r) — 7^ — r— tt^ — ~ — ; h 



+ J n+ i_F(/cro)Fn+i-F(A;r) 



Jl +F {kr Q ) + F„ 2 +F (A;ro) 
J n+ i_ F (A;ro)F n+ i_F(A;r) - F n+1 _^(A;ro) J n+1 - F (kr) 



Jn+i- F {kr ) + Y^ +1 _ F {kr ) 



E 

71=0 L 



Jl+ F {kr ) 



J 2 n+F (kr) + Y* +F (kr) 

J 2 n+F ( kr °)+ Y n +F ( kr °) 



n+l-F 



(kro) 



^ +1 _ F (fcr ) + i; 2 +1 _F(^o)J ' 

(15) 



In the absence of the magnetic flux in the tube we have 

S(Ax, fcr Q )|* =0 = S + Si (Ax, kr ), 

where 

oo 

Sq = J 2 (Ax) + 2^J 2 (Ax) = l, 

n=l 

and a correction term due to the finite thickness of an empty tube: 



(16) 
(17) 



Si (Ax, Ax ) 



T n w n \ Jo{kro)Y Q (kr) - Y (kr )J (kr t 



+ 2^ J n {kr )Y n {k 

n=l 

Jo(Ax 



J^kr ) + Y \kr ) 
J n (kr )Y n (kr) - Y n (kr )J n (kr) 



J*{kr Q ) + Y*{kr 



J*(kr ) + Y*(kr ) 
00 



2n \ J 2 (kr) + Y \kr) , o ^ ^ (Ax) + F n 2 (At) 



71=1 



J2(fcr ) + F n 2 (fcro) 



18) 



Thus, vacuum energy density (fTOi) depends on F fTT2|) . i.e. it is periodic in flux $ 
with a period equal to 2ne~ 1 . Moreover, relation (fTUj) is symmetric under substitution 
F — > 1 — F, vanishing at F — > (F — > 1) and, perhaps, attaining its maximal value at 
F = 1/20. Relations CHD and flE} are simplified at F = 1/2: 



2kr 



S (Ax)| 



<3>=7re _1 



sm r, 



r 



(19) 



1 At least, this is certainly true in the case of the singular vortex both for the Aharonov-Bohm [S] and 
the Casimir-Aharonov-Bohm [TU] - [T3] effects. 
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and 



Sx(kr, A;r )|$ =7re -i = 

~ JLi (kr ) \y* (At) - J 2 , (kr)] -2J n+k (kr )Y n+k (kr )J n+h (kr)Y n+h (kr) 



n=0 



J'Mk^+Y'ikro) 



(20) 



Since it is hardly possible to evaluate sums in (115]) and ( II 8ft analytically, our further 
analysis will employ numerical calculation. In the following we restrict ourselves to the 
case of F — 1/2, when the expression for the vacuum energy density takes form 



J dkk (k 2 + m 2 ) 1/2 G(kr, kr )- 
o 

oo 

-(£-1/4) A J dkk(k 2 + m 2 )~ 1/2 G(kr, kr, 
o 

G(kr,kr ) = S(kr, fcro)|$ =7re -i - S(kr, fcr )|$= . 



where 



(21) 



(22) 



3 Numerical evaluation of the vacuum energy density 

Following Ref. [T3j we rewrite ( 12 ip in the dimensionless form 

s / x ,„ i .s i , Q_(mro,mr) . , 

r 3 e ren = a + (mr Q , mr) - (£ - 1/4 r 3 A — ^ 23 

r 

where 



a±(mro,mr) = — dz z 
2ti 



z z + 



^rnr y 



±1/2 



G(z,Xz) 



(24) 



and A = r /r, A G [0, 1]. Let us point out some analytical properties of the integrand 
function in (T23|) : it vanishes at the edge of the vortex 

(25) 



lim G(z, Xz) = 0; 

A— >1 

at large distances from the vortex the case of a singular vortex is recovered 

lim G(z, Xz) = So(z)\$ =7re -i - S ; 

A->0 



at small values of z one gets 



G(z,Xz)\ z ^ = -[\n(X)/\n(Xz)} 2 . 



(26) 



(27) 



Numerical analysis indicates that in the calculation of function G(z, Xz) one can use 
series in ffT8]) and ([20]) with finite limits, namely for calculation G(z, Xz) at point z — z' 
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G(z,Xz) 



h=0. 



h=0. 




Figure 1: Behavior of G(z, Xz) at different values of A. 



it is enough to cut off the summation limits by some value n that can be found from 
condition 

< 5(A), 5(A) < 1(T 17 , (28) 

where N is a big number iV ~ 10 2 , n < N . It can be shown that the envelope of G(z, Xz) 
is exponentially decreasing function at large z, see Fig.l. So, for the finite-thickness 
magnetic vortex we can compute values of dimensionless quantity r 3 e ren ([2"3")) for different 
values of A. To do this, we have to be able to perform integration in f T2"3"j) with high 
precision. We make it in a following way. 

As one can see from Fig.2, the function G(z, Xz) is negative from z = to the first 
function root at z = Z\ {z\ ^ 0). So, the appropriate integral in ([211 is negative. The 
subsequent roots are denoted by z 2 , z 3 , etc. Because of decreasing character of the en- 
velope function the integral from z\ to z 3 will be positive. It is useful to define a period 
of function G(z, Xz) as an interval between two next to neighboring roots, i.e. from z\ to 
Zs, from z% to Z5, and so on. Then the full integral in (121]) will be a sum of the negative 
integral from z = to z = Z\ and a multitude of positive integrals over subsequent peri- 
ods. In the case of sufficiently small transverse size of the tube (mr <C 1) the integrals 
over some finite number of first periods may be negative but thereupon they become and 
remain positive also. 

For small z we make a direct integration of function G(z, Xz) over periods. For large 
z we make integration for each period separately. To do it we create a table of values 
of function G(z, Xz) for a separate period and replace this function by a more simple 
function in the form 

e -bz < z 2\ 

G int (z,Xz) = a - q sin(kz + jlnz + 4> ), (29) 
* ) 

where sine function ensures that roots of Gi nt (z, Xz) coincide with roots of G(z, Xz); A q (y) 
and B q (y) are g-degree polynomials, q can be 3,4 or 5; all unknown parameters can be 
found from an interpolation procedure. We allow a relative error of interpolation to be 

G int (z, Xz) - G(z, Xz) 
G(z,Xz) 



G(z,Xz)\ n - G(z,Xz)\ N 
G(z,Xz) 



< 10 



-8 



(30) 
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Figure 2: The location of roots of G(z, Xz) at A = 0.7. 



for each period. The function Gi n t(z, Xz) can be immediately integrated with the required 
accuracy. 

With the help of the above procedure we obtain a table of contributions from inte- 
gration over each period, extrapolate this table to infinity, and after that we find the 
full integral in fl2T]) as a sum of the negative integral over first period(s), a multitude of 
positive integrals over periods and an interpolation term. 

For a + function we estimate the relative error of the obtained result as 0.1%. It 
should be noted that nearly 95 % of the integral value is obtained by direct calculation 
and only nearly 5% is the contribution from the interpolation. The integration in a_ 
function is performed more quickly and with a higher accuracy, as compared to the case 
of a + function, because of its more rapid decreasing at large distances. In this case the 
contribution from the interpolation can be 10 _3 % from the final value of integration. 

Dimensionless quantity r 3 e ren fT2"3"j) can be interpreted as a function of two dimension- 
less parameters, mr and mr. Using the above described procedure, we calculate a + and 
ct_ functions at fixed values of mro (mro = 1; 10 _1 ; 10~ 2 ) at some set of points of dimen- 
sionless distance from the center of the vortex. This allows us to obtain coefficients^) of 
the interpolation function that is found in the form 



a±(xo, x) = [±e 



-2i x l=Fl/2 



X 



X 



X 



Pt{x 



Xq) 



X" 



Q 



X 



X ^> Xq, 



(31) 



where x = mr, xq = mro and P^, Qn — are polynomials of n-th order. First factor in the 
square bracket in fl3Tj) describes the large distance behavior of the appropriate functions 
in the case of the singular vortex [13] , second factor in the square bracket is an asymptotic 
at small distances from the edge of the tube, and the last factor is the intermediate part 
of the function. Since the flux tube is impenetrable, the a± functions are zero at x < x . 
Behavior of the dimensionless a± functions is presented on Fig.3. 
We define function 



a_(£o, x ) 



r 3 A 



a_(x , x) — x 



dct-(xo,x) 2 <9 2 a_(xo,a;) 



dx 



+ x 



dx 2 



(32) 



2 These coefficients arc different for different values of 



mr . 
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Figure 3: Behavior of the a+(xo,x) (positive) and the a-(xo,x) (negative) functions for 
the case of a) xo = 1, b) x = 10 _1 , c) xo = 1CT 2 . The variable x (x > xq) is along the 
abscissa axis. 



and present its behavior on Fig.4. 

Now we can construct the dimensionless vacuum energy density at different values of 
the coupling to the space-time curvature scalar: 



r 3 e T 



a + (x ,x) — (£ — l/4)a_(x , x). 



(33) 



Its behavior is presented on Fig.5 and Fig. 6. The case of the singular magnetic vortex 
[T3] is presented on Fig. 7. 

The analytical form of the vacuum energy density allows us to obtain the total induced 
vacuum energy 



2-7T OO 



E 



dip I e rPn rdr = 27rm 



(34) 



The integral over the a_ function can be taken by parts, yielding 



OO 

/ 



a_(xo,x) da_(x ,x) 
ax = — 



x z 



dx 



(35) 



X=XQ 



In this respect the question about the small distance behavior of the a_ function is very 
important. We have made a numerical calculations at small distances from the tube 
(x — xo ~ 10~ 6 £o) an d confirm the quadratic behavior near the edge of the tube (13T]) 
lim ac-(x ,x) ~ (x — x ) 2 . So, quantity (|35|) is zero, and the <S_ function affects only 



X— >XQ 



4-10 



2-10 



-2-10 




4-10 



-4-10 



-8-10 




2-1CT 



-410 



■10 




Figure 4: Behavior of the a_(x ,a;) function for the case of a) x = 1, b) x = 10 1 , c) 
xq = 10 -2 . The variable x (x > xq) is along the abscissa axis. 
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Figure 5: The dimensionless vacuum energy density r 3 e ren (xo, x) at different values of the 
coupling to the space-time curvature scalar for the case of a) x = 1, b) x = 10 _1 . The 
variable x (x > xq) is along the abscissa axis. 



the local properties of the vacuum energy density. The total vacuum energy is defined 
exclusively by the a + function and is independent of £: 



E = 2nm [ a+(Xo ' x) dx. (36) 
J x 2 

XO 

The total vacuum energy (in 2irm units) of the impenetrable flux tube is 6.94 • 1CT 10 , 
1.65 • 10 -4 , and 1.06 • 10 -2 for the case of xq — 1, xq — lO^ 1 , xq = 10~ 2 correspondingly. 
It should be noted that the total energy is infinite [13] in the case of a singular magnetic 
vortex. 

The induced vacuum energy density in the units of m 3 is presented on Fig. 8. The 
results for the case of £ = 1/4 are vanishingly small as compared to cases of other values 
of £, so they are not visible on Fig. 8. 




Figure 6: The dimensionless vacuum energy density r 3 e ren (xo, x) at different values of the 
coupling to the space-time curvature scalar for the case of Xq = 10~ 2 . The region in a 
rectangle on the left figure is seen in the scaled-up form on the right figure. The variable 
x (x > xq) is along the abscissa axis. 



9 




Figure 7: The dimensionless vacuum energy density r 3 e ren (xo, x) at different values of the 
coupling to the space-time curvature scalar for the case of the singular magnetic vortex. 

4 Discussion 

In the present paper we have considered the energy density which is induced in the vacuum 
outside a magnetic flux enclosed into an impenetrable tube of finite radius tq. Whereas 
the induced vacuum energy density is divergent at small distances as r -3 when the radius 
is neglected (r = 0), see Fig.7, it becomes finite when the radius is taken into account. 
A very characteristic feature is the appearance of oscillations in the vicinity of the tube, 
see Fig. 5 and Fig.6. Another peculiarity is that curves corresponding to different values 
of £ are symmetric with respect to the curve corresponding to £ = 1/4, the latter yielding 
the minimal absolute values of the vacuum energy density, see also Fig. 8. The maximal 
values of the vacuum energy density are becoming hardly observable for mr > 1, however, 
they are quite conspicuous for mro < 10~ 2 . This result was obtained earlier [TJ], but a 
completely new result concerns the behavior at large distances from the tube (up to 
200 r ), as well as at different values of £. It should be noted that the vacuum energy 
density in the vicinity of the tube is negative at £ < 1/4, including the important cases 
of conformal coupling £ = 1/8 (see (pQ) at d — 2) and minimal coupling £ = 0. 

Since the vacuum energy density is finite, the total vacuum energy, see (1511) . is finite 
as well. We show that the latter is positive and independent of £. Being negligible in 
the case of mro ^ 1 , it produces an appreciable effect of order of 10 _2 m in the case of 
mro — 10 -2 . 
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Figure 8: The vacuum energy density e ren (x ,x) (in m 3 units) at different values of the 
coupling to the space-time curvature scalar for the case of a) x = 1, b) x = 10 _1 , c) 
x = 10" 2 . 
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